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t""* '■ ABSTRACT 



We present the most sensitive 3 mm-survey to date of protoplanetary disks carried in the Taurus- Auriga star forming 
region (average rms of about 0.3 mjy), using the IRAM PdBI. With our high detection rate of 17/19, we provide the 
first detections at wavelengths longer than about 1 mm for 12 sources. This enables us to study statistically the mm 
SED slopes and dust properties of faint disks and compare them to brighter disks using a uniform analysis method. 
With these new data and literature measurements at sub- millimeter and millimeter wavelengths, we analyze the dust 
properties of a sample of 21 isolated disks around T Tauri stars in the Taurus-Auriga star forming region. Together 
with the information about the disks spatial extension from sub/mm-mm interferometric studies, we derive from the 
observed sub-mm/mm spectral energy distribution constraints on the dust opacity law at these wavelengths, using 
two-layer flared disk models and a self-consistent dust model that takes properly into account the variation of the dust 
opacity with grain growth. We find evidence for the presence in the disk midplane of dust particles that have grown to 
sizes as large as at least 1 millimeter in all the disks of our sample, confirming what was previously observed on smaller 
brighter objects. This indicates that the dust coagulation from ISM dust to mm-sized grains is a very fast process in 
protoplanetary disks, that appears to occur before a young stellar object enters the Class II evolutionary stage. Also, 
the amount of these large grains in the disk outer regions is stationary throughout all the Class II evolutionary stage, 
indicating that mechanisms slowing down the dust inward migration are playing an important role in the Taurus- Auriga 
protoplanetary disks. Another result is that the spectral index between 1 and 3 mm for the 6 faintest disks in our sample 
is on average smaller than for the brighter disks, indicating either that these fainter, yet unmapped, disks are spatially 
much less extended than the brighter spatially resolved disks, or that fainter disks have typically larger dust grains in 
their outer regions. Considering that these fainter disks are more representative of the bulk of the disk population than 
£N| ■ the brighter ones, this may have important consequences for the theories of planetesimal formation and disk formation 

and evolution. Finally, we investigate the relations between the derived dust properties, namely dust mass and grain 
^\ ] growth, and the properties of the central star, like its mass, age and mass accretion rate. 
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1. Introduction ferent wavelengths. Throop et al. (|2001D observed in optical 

and near- infrared the translucent edge of the 114-426 pro- 
Circumstellar disks play a fundamental role in the physical toplanetary disk in the Orion Nebula Cluster, and derived 
processes involved in star and planet formation. In these f rom the extinction curve of the background nebular emis- 
systems part of the material (gas and dust) loses angu- s i on a i ower limit of a few microns for the dust grains in the 
lar momentum and accretes onto the central forming star, very oute r regions of the disk, about 500 AU away from the 
while another part of it may give birth to a planetary sys- cen tral star. Another method to probe the presence of large 
tern. Although the content of dust is only a small fraction grains m t ^ e disk is through the shape and intensity of the 
of the overall material in a circumstellar disk, dust grains iq. ancl 20-^m silicate features. The data indicate a large 
are crucial elements in the early stages of planet forma- var i e t y of silicate profiles, ranging from strongly peaked 
tion. According to the core-accretion scenario (Safronov & silicate bands and steeply rising spectral energy distribu- 
Zvjagina H9691 Pollack et al. |1996p, growth from an ini- t ions (SEDs) to " boxy" silicate profiles and flatter SEDs 
tial ISM population of submicrometer-sized dust grains to (iso, Malfai t et a l. gMj Spitzer Space Telescope, Kessler- 
planetesimals of 1 - 100 km sizes is believed to be the key Silacci et al. CHU). The boxy features with low feature-to- 
mechanism of early planet formation. These planetesimals continuum ratios are interpreted as grain growth to micron 
"embryos" then continue to grow leading to the formation s i ze (Bouwman et al. I2001|) . 
of at least terrestrial planets and possibly the cores of giant 

planets. However, the techniques outlined so far can only probe 

Observational evidence of grain growth from the ISM dust grains as large as some /an in the disk surface layers, 

dust has been obtained from a variety of techniques at dif- They are not sensitive to larger grains in the disk mid- 
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plane, the region in which planet formation is supposed to 
occur. In order to probe larger grains in the midplane, ob- 
servations at millimeter wavelengths are needed. Beckwith 
& Sargent (|199ip showed that T Tauri stars have shallow 
SEDs at sub- millimeter wavelengths. Under the assumption 
of optically thin emission at these frequencies, this implies a 
dust opacity dependence on wavelength (K mm oc A - ' 5 ) much 
flatter than in the ISM (/3ism = 1-7), which is naturally in- 
terpreted in terms of grain growth (see e.g. Draine 2006). 
However, this interpretation of the disk SEDs is not unique, 
since the same data can be explained by very small opti- 
cally thick disks. To break this degeneracy and sort out 
the effect of potentially large optical depth it is neces- 
sary to spatially resolve the disks to determine their actual 
sizes. Furthermore, in this context observations at millime- 
ter wavelengths are very useful, since at these lower fre- 
quencies the impact of optically thick disk inner regions 
to the total emission is expected to be lower. Therefore, 
to actually probe the dust properties in protoplanetary 
disks, one needs to combine the determination of the sub- 
millimeter /millimeter SED with information on the disk 
extension from high-angular resolution interferometric ob- 
servations (see e.g. Testi et al. [2001). 

In the last years, several sub-mm/mm interferometric 
observations have been carried out to investigate dust grain 
growth in protoplanetary disks. Wilner et al. (2000) re- 
solved the disk around the TW Hya pre-main-sequence 
(PMS) star at 7 mm using the Very Large Array (VLA). 
Extensive modelling of the SED of this source has shown 
that the dust grains in the outer parts of the TW Hya 
disk have grown to at least ~ 1 cm (Calvet et al. 2002). 
An analogous result has been obtained for the disk around 
CQ Tau (Testi et al. 12003]) . More recently, Rodmann et 
al. (|2006j) resolved 10 disks in the Taurus- Auriga star form- 
ing region and found clear evidence of grain growth in 4 of 
them. Lommen et al. (2007) resolved 1 disk in Chamaeleon 
and 4 in Lupus with the Submillimeter Array (SMA) at 
1.4 mm and with the Australia Telescope Compact Array 
(ATCA) at 3.3 mm, and found clear evidence of dust grain 
growth to sizes of a few millimeter for 4 of them. Finally 
Schaefer et al. (|2009[) observed a sample of 23 low mass 
PMS stars (with spectral types of K7 and later) in Taurus- 
Auriga at 1.3 mm and 2.6 mm with the Plateau de Burc 
Interferometer (PdBI); they detected only 8 sources at 1.3 
mm and 6 at 2.6 mm, and found evidence of grain growth 
for the 3 disks that they could spatially resolve. 

In this paper we present our analysis on 21 protoplane- 
tary disks around T Tauri stars in the Taurus- Auriga star 
forming region without stellar companions in the range 
0.05" — 3.5". For 11 of these objects, mainly faint disks with 

3mm < 100 mjy, we have obtained new data at ~ 3 mm 
with the Plateau de Bure InterferometeiQ (hereafter PdBI) 
with an average rms of about 0.3 mjy. Together with the 
data already present in the literature at sub/mm- and mm- 
wavelengths, and with the information obtained in the last 
years on the disks spatial extension from high-angular reso- 
lution interferometric observations at mm wavelengths, we 
investigate the dust properties in the disks, namely grain 
growth and dust mass, and their relation with the proper- 
ties of the central star. 



1 The Plateau de Bure Interferometer at IRAM is supported 
by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 



In Section [5] we present our new PdBI data, the proper- 
ties of our sample, the method we used to estimate the main 
stellar physical quantities, and the data we used for our 
analysis. In Section [3] we describe the disk models adopted 
for the analysis of the disks sub-mm/mm SED. In Section 
[4] we show and discuss the results of our study in terms of 
dust grain growth and dust mass, whereas in Section [5] we 
summarize our main findings. 



2. New 3mm observations and sample properties 

2.1. New PdBI observations 

We observed a total list of 19 targets; 16 of them were se- 
lected to be relatively faint, with 15 mJy< Fi. 3mm < 100 
mjy and with no known detections beyond about 1mm; 
the other 3 (HL, DO and DR Tau) are brighter sources 
that were observed for references purposes. The observa- 
tions were carried out with PdBI between July and August, 
2007. Owing to antenna maintenance work, all observations 
were carried out in a subarray of the six-element interfer- 
ometer. The antennas were arranged in very compact con- 
figurations that provided sensitive baselines from 15 m to 
80 m. The dual-polarization receivers, which were observ- 
ing in the 3mm band, had typical receiver temperatures of 
about 35 K and were providing image sideband rejection 
values better than 10 dB. To achieve maximum sensitivity 
in the continuum, the spectral correlator was adjusted to 
cover a total effective bandwidth of about 2 GHz. 

The observing procedures were set up to fill in gaps 
in the scheduling without any specific tuning wavelength 
within the 3mm band. Visibilities were obtained using on- 
source integration times of 22 min interspersed with 2 min 
calibrations on a nearby calibrator. Each star was observed 
on one or more occasions for a minimum of 40 min and a 
maximum of 2 hr on-source. The atmospheric phase stabil- 
ity on the baselines was always better than 40° , consistent 
with seeing conditions of l"-3", typical for summer. The ab- 
solute flux density scale was calibrated on 3C84, 0528+134, 
MWC349 and Mars, and was found to be accurate to 10%. 
The receiver passband shape was determined with excellent 
precision on a strong calibrator like 3C84. 

We used the GILDAS package for the data reduction 
and analysis. The data were calibrated in the antenna based 
manner. The continuum visibilities were gridded with nat- 
ural weighting and no taper to maximize the sensitivity. 
Since none of the target stars was found to be resolved in 
synthesized beams of 2"-4", a point source model was fit- 
ted to the calibrated visibilities of each star to estimate the 
flux density of the continuum in the 3mm band. The results 
are summarized in Table [T] Among the 19 targets, 17 were 
detected, including 14 of the 16 faint sources, thanks to our 
good sensitivity; 12 of these are first detections longward of 
about 1mm, strongly expanding the available database of 
3mm data towards fainter disks. The 2.6mm detections of 
LkHa 358 and GO Tau were recently published by Schaefer 
et al. (|2009)) . 

In Sections 12. 2\ 12. 31 and 12. 41 we describe the final sample 
used for our analysis of dust properties, the adopted method 
to estimate the stellar properties and mass accretion rates 
summarized in Table [21 whereas in 12.51 we report the sub- 
mm/mm data that we collected from the literature for our 
analysis. 
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Object a (J2000) 8 (J2000) v X F v rms Comments 

(GHz) (mm) (mjy) (mjy) 



CW Tau 


04:14:17.0 


28: 


:10 


:56 


.51 


84.2 


3.57 


3.44 


0. 


.31 




CX Tau 


04:14:47.8 


26: 


:48 


:11. 


.16 


86.2 


3.49 


1.01 


0. 


.13 




DE Tau 


04:21:55.6 


27: 


:55 


:05. 


.55 


92.6 


3.23 


3.32 


0. 


.25 




DK Tau 


04:30:44.3 


26: 


:01 


:23. 


.96 


84.2 


3.57 


4.70 


0. 


.20 


2.5" binary 


DO Tau 


04:38:28.6 


26: 


:10 


:49. 


.76 


84.2 


3.57 


13.80 


0. 


.29 




DP Tau 


04:42:37.7 


25: 


:15 


:37. 


.56 


84.2 


3.57 


< 0.78 


0. 


.26 


no sub-mm info 


DR Tau 


04:47:06.2 


16: 


:58 


:43. 


.05 


100.8 


2.97 


13.90 


0. 


.47 




DS Tau 


04:47:48.6 


29: 


:25 


:10. 


.96 


100.8 


2.97 


2.96 


0. 


.21 




FM Tau 


04:14:13.5 


28: 


:12 


:48. 


.84 


100.8 


2.97 


2.65 


0. 


.24 




FP Tau 


04:14:47.4 


26: 


:46 


:26. 


.65 


100.8 


2.97 


11.80 


0. 


.26 


no sub-mm info 


FY Tau 


04:32:30.5 


24: 


:19 


:56. 


.69 


100.8 


2.97 


1.90 


0. 


.28 


no sub-mm info 


FZ Tau 


04:32:31.8 


24: 


:20 


:02. 


.53 


100.8 


2.97 


2.89 


0. 


.24 




GK Tau 


04:33:34.6 


24: 


:21 


:06. 


.07 


100.8 


2.97 


< 0.96 


0. 


.32 


2.5" binary 


GO Tau 


04:43:03.1 


25: 


:20 


■17. 


.38 


86.2 


3.49 


4.05 


0. 


.34 




Haro 6-28 


04:35:56.9 


22 


M 


:36 


.63 


87.0 


3.45 


0.92 


0. 


.18 


0.7" binary 


HL Tau 


04:31:38.4 


18: 


:13 


:57. 


.37 


84.5 


3.53 


43.40 


0. 


.68 


flat spectrum 


HO Tau 


04:35:20.2 


22: 


:32 


:13. 


.98 


87.0 


3.45 


2.03 


0. 


.53 




LkHa 358 


04:31:36.2 


18: 


:13 


:43. 


.20 


84.5 


3.53 


2.46 


0. 


.22 


no sub-mm info 


SU Aur 


04:55:59.4 


30 


:34 


:10. 


.39 


84.5 


3.53 


2.88 


0. 


.19 





1) Reason why the source has not been considered in the analysis (see text for more details). 



2.2. Sample 



Our final sample of 21 sources consists of all young stellar 
objects (YSO) in Taurus- Auriga catalogued in Andrews & 
Williams (2005) that fulfill the following criteria: (1) a Class 
II infrared SED, in order to avoid contamination of sub-mm 
fluxes by a residual envelope; (2) information on the central 
star through optical-NIR spectroscopic/photometric data, 
necessary to calculate self-consistent disk SED models; (3) 
one detection at ~ 3mm (either from our new PdBI obser- 
vations or from the literature) and at least one detection in 
the 0.45 mm < A < 0.85 mm spectral region. This selection 
criterion was chosen in order to probe a broad enough spec- 
tral window to efficiently constrain the optically thin part 
of the sub-mm/mm spectral energy distribution (SED); (4) 
no evidence of stellar companions in the 0.05 — 3.5" range in 
angular separation. This range corresponds to about 5 — 500 
AU in projected physical separation at the Taurus- Auriga 
star forming region distance, here assumed to be 140 pc for 
all the sources of our sample (Bertout et al. U999j) . The limit 
of 3.5 arcsec ensures that mm interferometric flux measure- 
ments, of typical resolution 1 arcsec, are not contaminated 
by the companion's disk, while binaries closer than 5 AU 
should have outer disk properties similar to those of single 
stars on the scales probed by (sub-)mm data. 

Only 21 sources catalogued by Andrews & 
Williams (|2005|) fulfill the above 4 criteria. They are 



listed in Table [H and include 11 sources from our PdBI 
sampltQ. 

Figure 1 reports histograms that describe the com- 
pleteness level of our final sample with respect to all 
Taurus- Auriga Class II YSOs catalogued in Andrews & 
Williams (2005), and fulfilling the "isolation" criterion (4). 
Our sample contains 63% of these "isolated" class II show- 
ing a flux at ~ 0.85mm greater than 100 mjy, and 31% of 
the sources with lower 0.85mm fluxes. Therefore, our sam- 
ple is not complete even for the brightest YSOs at 0.85mm. 
This is due to a lack of observations at ~ 3mm for these 
sources. However, thanks to the high sensitivity of our new 
PdBI observations (see Section |2~T]) our sample enables us 
to study statistically for the first time the dust properties 



2 The 8 sources dropped from the PdBI sample are HL Tau 
(a flat spectrum source that retains an envelope; Padgett et 
al. [T999|) : DP Tau, FP Tau, FY Tau and LkHa 358 that have 
no available submm detections for SED fitting; and DK Tau, 
GK Tau and Haro 6-28 which are binary systems with an- 
gular separations of 2.5" (Simon et al. I1992|l . 2.5" (Reipurth 
& Zinnecker ll"993j) and 0.7" (Leinert et al. |1993|) respectively. 
We keep in our sample the wide binaries DS Tau and HO Tau 
(angular separations of 7.1" and 6.9" respectively, Mathieu et 
al. 11994 [) ; they do not show any emission from the compan- 
ions from our new 3 mm interferometric observations. The other 
dropped multiple systems from the AW05 catalog are DH Tau 
(2.3" binary, Itoh et al. l2005[) , GG Tau (quadruple system with 
a circumbinary ring of inner and outer radii of 180 and 260 AU 
respectively, Pinte et al. [20071), and UY Aur (0.9" binary, Leinert 
et al. ri993|) 
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of faint disks (i.e. with Fi. 3mm < 100 mJy, F , S5mm < 100 
mjy) and compare them to the brighter ones using a uni- 
form analysis method. In terms of the stellar properties, it is 
important to note that the Andrews & Williams (2005) cat- 
alogue includes only a few very low mass PMS stars: only 
stars with an estimated mass higher than about 0.2 M Q 
have been detected because of the sensitivity limits of the 
current facilities at sub-mm and mm wavelengths that make 
the detection of disks around very low mass YSOs very diffi- 
cult. The sensitivity limitation is even more severe at 3mm, 
therefore our final sample includes 58% of all the class II 
isolated sample of Andrews & Williams (2005) with an esti- 
mated mass greater than 0.4 Mq, but only two sources with 
M* < 0.4 M Q (out of 10 in Andrews & Williams HMD . 

2. 3. Stellar properties 

Spectral types for all the sources in our sample are taken 
from the literature (Kenyon & Hartmann 119951 Briceno et 
al. [TMSj) and cover a range from G2 (SU Aur) to M2.5 
(CX Tau). Spectral types of M0 and earlier have been con- 
verted to effective temperatures with the dwarf tempera- 
ture scale of Schmidt-Kaler (|1982p . whereas for types later 
than M0 we adopted the temperature scale developed by 
Luhman 119991 Typical uncertainties for these sources are 
~ 1 — 1.5 in spectral sub- types or ~ 100 — 150 K in tem- 
perature. 

The stellar luminosities of all the pre-main-sequence 
(PMS) stars were computed from their J-band flux (2MASS 
All Sky Catalog of point sources, Cutri et al. 2003} to min- 
imize contamination from UV and IR excess emission. For 
the bolometric corrections we adopted the dwarf values 
from Kenyon & Hartmann (|1995[) that are considered to 
be satisfactory approximations for these young sources in 
the J-band (Luhman UM§§ ■ 

In order to estimate the amount of extinction toward 
these young sources, we calculated the color excesses using 
intrinsic colors provided by Kenyon & Hartmann (j 1995|) for 
G2-K7 spectral types and by Leggett (fT9l?2)) for M0-M2.5 
types. To ensure that the color excess reflects only the ef- 
fect of reddening, minimizing the emission due to accretion, 
the typical selected colors are between the R and H bands. 
Since for some of the objects at later spectral types in our 
sample R — I measurements are not available, we dered- 
dened the J - H and H - K s colors from the 2MASS All 
Sky Catalog of point source^] to the locus observed for clas- 
sical T Tauri stars (CTTS) by Meyer et al (fT997j) . following 
the method described in Briceno et al. (2002). Extinctions 
arc finally calculated adopting the extinction law of Rieke 
& Lebofsky (PSS]) . 

The computed photospheric luminosities cover a range 
from 0.16 L Q (HO Tau) to 10.6 L Q (SU Aur). Considering 
the uncertainties in the photometry, reddenings, bolometric 
corrections and distance, the typical errors in the bolomet- 
ric luminosities are ±0.08 — 0.13 in logi*. Slight differ- 
ences from the values obtained by other authors are typi- 
cally within the uncertainties and are not significant for the 
purposes of this paper. 



For this analysis the JHK S magnitudes have been 
transformed from the 2MASS photometric system to the 
Johnson-Glass one using the color transformations reported in 
Carpenter (|200Tj) and Bessell & Brett (fl988|) . 



10.00 r--- 



1.00 



0.10 



0.01 




6000 



5000 



4000 



3000 



(K) 



Fig. 2. H-R diagram for the sources of our sample. The 
dashed lines correspond to the isochrones for ages in Myr 
as labeled at the right end of the lines from the [T999 PMS 
evolutionary models, while the solid ones represent the evo- 
lutionary tracks from the same models for PMS stars with 
masses as labeled at the top of the evolutionary track lines. 
In this diagram the evolutionary tracks start from an age 
of 0.1 Myr. The errorbars for the (T e s, L+) values of our 
objects are not shown, the uncertainties in the photomet- 
ric parameters are typically ~ 0.08 — 0.13 in logL* and 
~ 100 — 150 K in temperature. 



Given the effective temperatures and photospheric lu- 
minosities, we placed the PMS stars on the H-R diagram 
and derived stellar masses and ages adopting the theoretical 
tracks and isocrones of Palla & Stahler Q1999P (see Figure 
[5]) . One important source of uncertainty for stellar masses 
and ages is given by the spread of values obtainable by 
using different PMS evolutionary models. For example, if 
compared to the ones obtained with the Baraffe et al. ( 1998) 
models, our adopted values of stellar masses and ages are 
typically lower by a factor of ~ 1.5 and ~ 2 respectively. 
Our choice of using the Palla & Stahler (|1999[) models is 
uniquely due to the more complete coverage of the HR di- 
agram plane, that allows us to use the same evolutionary 
models for all the sources of our sample. One should al- 
ways bear in mind the high uncertainty associated to these 
quantities, expecially to the stellar age (see Section [4*7T1 for a 
more deailed discussion about the estimates of YSO ages). 
According to the Palla & Stahler (1999) models the ranges 
spanned by our sample go from about 0.3 M Q (CX Tau) 
to 2.2 Mq (RY Tau) in stellar masses and from about 0.1 
Myr (UZ Tau E) to 17 Myr (DS Tau) in stellar ages. 

The main stellar quantities described here are summa- 
rized in Tabled 

2.4. Mass accretion rates 

Given the stellar luminosities and effective temperatures 
we calculated the stellar radii using the Stefan-Boltmann 
law. Once the stellar mass M* and radius i?* are set, the 
mass accretion rate M acc can be directly obtained from an 
estimate of the accretion luminosity L acc by 



= 1.25 



GM* 



(1) 
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Fig. 1. Histograms higlighting some properties of our selected sample. In all the histograms our sample is represented 
by black +grey columns, the black columns represent the disks in our sample for which we obtained new PdBI data at ~ 
3 mm, while the total columns (black+grey+white) include all the class II YSOs from the Andrews & Williams ( 2005) 
catalogue with no evidence of stellar companions in the 0.05 — 3.5" interval in angular separation, from which our sample 
has been selected (see text). Upper left histogram represents the distribution of the ^0.85 mm fluxes, including upper 
limits, for all the sources but two, 04113+2758 and GM Aur, that have not been observed at ^0.85 mm. Upper right and 
bottom histograms show the distribution of stellar spectral type and estimated stellar mass (see Section |2~3| respectively. 
Four class II YSOs are not included in these histograms because for them an estimate of the stellar spectral type is not 
available. These are 04301+2608, FT Tau, Haro 6-39 and IC 2087. 



where G is the gravitational constant and the factor of ~ 
1.25 is calculated assuming a disk truncation radius of ~ 
5 (Gullbring et al. 1X5551) . 

The estimates of L acc were obtained from spectroscopic 
detections of excess Balmer continuum emission in the lit- 
erature (see references in Table [2]). These estimates are 
obtained primarily from measuring the excess flux in the 
Balmer continuum shortward of 3646 A. The luminos- 
ity is then obtained by applying a bolometric correction 
calculated from shock models (Calvet & Gullbring 1998, 
Gullbring et al. [2000) or simplistic plane-parallel slabs 
(Valenti et al.QMS Gullbring et al. fTMg]) . 

Two sources in our sample, FZ Tau and GO Tau, 
do not have literature estimates of accretion luminosity 
from the excess Balmer continuum emission. We obtained 
the accretion luminosity for these two sources from low- 
resolution spectra from 3200-8700 A obtained with the 
Double Spectrograph at the 5m Hale Telescope at Palomar. 
The accretion luminosity was then calculated from the ex- 
cess Balmer continuum emission, following the method de- 
scribed in Herczeg & Hillenbrand (2008). 



2.5. Disks sub-mm and mm data from the literature 

In order to probe the sub-millimetre and millimeter SED 
we collected data of the dust continuum emission between 
~ 0.450 to ~ 7 mm from several works in the literature, 
listed in Table [3J These data, except for 7 mm measure- 
ments, have been used together with the new PdBI ob- 
servations at ~ 3 mm to constrain the disk properties by 
fitting the sub-mm/mm SED with the models described 
in Section [3l The 7-mm data is excluded because free-free 
emission from ionized gas contaminates disk emission, with 
free-free emission typically contributing ~ 20% of the 7mm 
flux (Rodmann et al. 120061) . Simultaneous cm observations 
are required to constrain the free-free emission at 7mm, 
which is variable on short timescales (< a few days, see 
Lommen et al. [2009). Only 2 disks in our sample (RY 
Tau and UZ Tau E) have been observed nearly simultane- 
ously at centimeter and millimeter wavelengths (Rodmann 
et al. I2006j) . without which the 7mm fluxes should be con- 
sidered only as upper limits for the dust emission. We there- 
fore did not include the 7mm data in our analysis, although 
we included the available 7mm data from the literature in 
Figure S] and verified a posteriori that the fluxes of our disk 
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Object 


ST 


T c ff 


Lace Ref. (I/acc) a 


Ri, Log M acc 




Age 






(K) (L Q ) 




(Jig) (M yr- 1 ) 


(M @ ) 


(Myr) 



AA Tau 


K7 


4060 


0.80 


0.03 


1 


1.8 


-8.56 


0.8 


2.4 


CI Tau 


K7 


4060 


1.01 


0.11 


2 


2.0 


-7.95 


0.8 


1.7 


CW Tau 


K3 


4730 


1.32 


0.05 


2 


1.7 


-8.58 


1.2 


6.6 


CX Tau 


M2.5 


3488 


0.48 


0.03 


2 


1.9 


-8.18 


0.3 


0.8 


CY Tau 


Ml 


3705 


0.55 


0.04 


1 


1.8 


-8.24 


0.5 


1.4 


DE Tau 


M2 


3560 


1.01 


0.07 


1 


2.6 


-7.69 


0.4 


0.2 


DL Tau 


K7 


4060 


0.89 


0.23 


2 


1.9 


-7.66 


0.8 


2.2 


DM Tau 


Ml 


3705 


0.24 


0.10 


2 


1.2 


-7.99 


0.5 


3.6 


DN Tau 


MO 


3850 


0.91 


0.02 


1 


2.2 


-8.55 


0.6 


1.1 


DO Tau 


MO 


3850 


0.79 


0.60 


1 


2.0 


-7.10 


0.6 


1.3 


DR Tau 


K7 


4060 


1.09 


0.72 


2 


2.1 


-7.12 


0.8 


1.5 


DS Tau 


K2 


4900 


0.90 


0.21 


1 


1.3 


-8.00 


1.1 


17.0 


FM Tau 


MO 


3850 


0.26 


0.10 


2 


1.2 


-8.13 


0.6 


6.3 


FZ Tau 


MO 


3850 


1.06 


0.46 


3 


2.3 


-7.14 


0.6 


0.8 


GM Aur 


K3 


4730 


1.26 


0.07 


1 


1.7 


-8.44 


1.2 


7.2 


GO Tau 


MO 


3850 


0.27 


0.01 


3 


1.2 


-9.14 


0.6 


4.8 


HO Tau 


M0.5 


3778 


0.16 


0.01 


2 


0.9 


-9.13 


0.5 


8.7 


IQ Tau 


M0.5 


3778 


0.73 


0.04 


2 


2.0 


-8.24 


0.6 


1.2 


RY Tau 


Kl 


5080 


6.59 


1.60 


4 


3.3 


-7.03 


2.2 


1.1 


SU Aur 


G2 


5860 


10.6 


0.55 


4 


3.2 


-7.52 


1.7 


2.2 


UZ Tau E 


Ml 


3705 


2.27 


0.31 


2 


3.7 


-7.02 


0.5 


0.1 



a) References. (1) Gullbring et al. (1998); (2) Valenti et al. () 1993(1 . these accretion luminosities have been multiplied by a factor 
(140 pc/160 pc) 2 to account for the current estimate of 140 pc for the distance to the Taurus-Auriga star forming region; (3) 
Herczeg et al. (2009); (4) Calvet et al. Ip004)) . 



models for all sources are equal or less than the measured 
flux at 7mm. Simultaneous 7mm and cm-wavelength obser- 
vations with the EVLA will extend the coverage of the disk 
SED to longer wavelengths, in order to better constrain the 
millimeter-wavelengths dust opacity and to provide infor- 
mation on dust grains larger than ~ 1 cm. 



3. Disk models 

We compare the observational sub-mm/mm data with the 
two-layer models of flared disks (i.e. in hydrostatic equi- 
librium) heated by the stellar radiation as developed by 
Dullcmond et al. (2001), following the schematization of 
Chiang & Goldreich (1997). These models have been used 
in the analysis of CQ Tau by Testi et al. (2003) and of a 
sample of 9 Herbig Ae stars by Natta et al. (2004j) and we 
refer to these papers for a more detailed description. 

To completely characterize a model of the disk, we need 
to specify the stellar properties (i.e. the luminosity L*, ef- 
fective temperature T e s, mass M* and distance d, assumed 
to be 140 pc for all the sources in our sample), some char- 
acteristics of the disk structure, namely the inner and outer 
radius (R ul and R ou t, respectively), the parameters Si and 
p defining the radial profile of the dust surface density as- 



sumed to be a power-la-wQ Sdust = £dust,i(-R/l AU)~ P , the 
disk inclination angle i (90° for an edge-on disk) and, lastly, 
the dust opacity. The millimeter SED is almost completely 
insensitive to some of the disk parameters, which cannot 
thus be constrained by our analysis. For example, the inner 
radius of the disk, i?i n , affects the emission only at near 
and mid-infrared wavelengths. 

The outer radius i? ou t is a parameter that can be probed 
by high-angular resolution interferometric observations. In 
our sample 13 disks have been spatially resolved, yielding an 
estimate for i? ut- It is important to note that different in- 
terferometric observations performed on the same disk can 
lead to very different values of i? Q ut (> 100 AU) because 



4 As shown by Hughes et al. (f2008|l and Isella et al. (2009), 
a tapered exponential edge in the surface density profile, physi- 
cally motivated by the viscous evolution of the disk, is able to ex- 
plain the apparent discrepancy between gas and dust outer radii 
derived from millimeter observations of protoplanetary disks. 
Here we want to note that when fitting the sub/mm- mm SED 
with an exponential tail instead of a truncated power-law, the 
dust opacity spectral index /3 (defined in Section 13. ip is un- 
changed, and disk masses vary by less than a factor of 2, even 
if the outer disk radius becomes much bigger (see, e.g., Table 
4 in Kitamura et al. I2002f> . So the simpler truncated power-law 
model is a justified approximation for the purposes of this paper. 
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A (mm) 


Source References" 


0.450 


x 


0.600 


9 3 


0.624 


4 


0.769 


4 


0.800 


9 3 


0.850 


1 5 


0.880 


5 


1.056 


4 


1.100 


2 3 


1.200 


(; 


1.300 


1, 7, 8, 9, 10 


1.330 


5 


2.000 


2, 11 


2.700 


7, 8, 12, 13 


3.100 


8 


3.400 


9, 14 


7.000 


15 



a) Refer ences. (1) Andrews & Williams ([2005]l; (2) Mannings & Emerson (T994} ; (3) Adams e t al. fl l990); (4) Beckwit h fc 
Sargent (TT99T]) ; (5) A ndrews & Williams (2007); (6) Althenoff et al. ([1994]) : (7) Dutrey et al. (fT996) ; (8) Jensen et al. (1996); (9) 
Koerner et al. (flM5jl; (1 0) Isella et al. (|2009]l : ( 11) K itamura et al. (|2002|l : (12) Mundy et al. {1996); (13) Schaefer et al. l[20T)9j> : 
(14) Ohashi et al. CSSSJl; (15) Rodmann et al. (12006)1 . 



of differences in the observations angular resolution, sen- 
sitivity, and methods for deriving the radius. Given these 
uncertainties in the determination of i? ut, we consider a 
range of possible values for i? ou t rather than a single-value 
estimate, as listed in Table [4] For the 8 objects in our sam- 
ple that have not yet been mapped with angular resolu- 
tions smaller than a few arcseconds which may potentially 
resolve the disk, a fiducial interval of 100 — 300 AU has 
been assumed^. This adopted range for R out comprises the 
vast majority of the disk outer radii distribution function 
as derived by Andrews & Williams (2007) through ~ 1" an- 
gular resolution sub-millimeter observations of 24 disks in 
Taurus- Auriga and p— Ophiucus. We will discuss the impact 
of the outer disk uncertainty on the results of our analysis 
in Section |H Here we want to note that with these values of 
i?out the disk emission at sub-mm/mm wavelengths is dom- 
inated by the outer disk regions (R > 30 AU) which are op- 
tically thin to their own radiation. An effect of this is that 
the millimeter SED does not depend on the disk inclina- 
tion angle i (for i <80°), thus decreasing the number of the 
model parameters used to fit the observational SED. This is 
certainly true for the spatially resolved disks for which we 
have an estimate for i? ou t- The 8 unmapped objects could 

5 We will refer to these sources (CW, CX, DE, DS, FM, FZ, 
HO Tau, and SU Aur) as "unmapped" for the rest of the paper. 
Note that they are all faint with F\ mm < 120 mjy. This probably 
explains why for these objects no observations at high-angular 
resolution have been attempted yet: the sensitivity of the current 
(sub-)mm interferometers would be not high enough to detect 
the low surface brightness of the outer disk regions. 



have optically thick inner disks that contribute significantly 
to the sub-mm/mm emission. Since the disk outer radius 
cannot be constrained by the SED alone, interferometric 
imaging that spatially resolves these disks is the only way 
to get information on the exact value of i? ut- Nevertheless, 
since the disks that have been spatially resolved so far show 
outer radii typically in the range 100 — 300 AU, we consider 
highly unlikely that a significant fraction of the unmapped 
disks have i? ut ^ 50 AU, which would make the disk opti- 
cally thick even at mm wavelengths. 

3.1. Dust opacity 

The last quantity that defines the disk model is the dust 
opacity, which depends on grain sizes, chemical composi- 
tions, and shapes (e.g. Miyake & Nakagawa 119931 Pollack 
et al. 119941 Draine et al. [2006). In order to extract from 
the sub-mm/mm SED quantitative estimates for the dust 
grain sizes and the dust mass in the disks, assumptions on 
the chemical composition and shape of the dust grains have 
to be made. 

For our discussion we considered porous composite 
spherical grains made of astronomical silicates, carbona- 
ceous materials and water ice (optical constants for the 
individual components from Weingartner & Draine 2001, 
Zubko et al. 119961 Warren [ 19841 respectively) adopting frac- 
tional abundances used by Pollack et al. (19911) ( see caption 
of Figure [3]). We considered a population of grains with a 
power-law size distribution n{a) oc a~ q between a mini- 
mum and a maximum size, a m i n and a max , respectively. In 
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the ISM, a m ; n is a few tens of A, a max ~ 0.1 — 0.2 /im, 
and q = 3.5 (Matins et al. [WT1 Draine & Lee [T9"8"4l by 
analysing extinction and scattering of starlinght from the 
interstellar dust). As discussed e.g. in Natta et al. (2004) 
in a protoplanetary disk, because of grain processing, one 
expects much larger values of a max than in the ISM, and 
a variety of possible q values, depending on the processes 
of grain coagulation and fragmentation that occur in the 
disk. Experimental studies of fragmentation from a single 
target have found values of q ranging from w 1.9 for low- 
velocity collisions to « 4 for catastrophic impacts (Davies 
& Ryan 1990), but the power-law index for fragments from 
shattering events integrated over a range of target masses 
need not be necessarily the same as in the case of a single 
target. In the disk, if the coagulation processes from the 
initial ISM-grains population dominate over the fragmen- 
tation ones, then a value for q that is lower than the ISM 
one is expected (the opposite should occur in the case of 
fragmentation playing the major role). 

In the adopted two-layer disk model, the disk is divided 
into the surface and midplane regions, and therefore two 
different opacity laws are considered. Because of the as- 
sumed vertical hydrodynamical equilibrium and dust set- 
tling, in a protoplanetary disk the content of mass in the 
midplane is much higher than in the surface. At millimeter 
wavelengths, where the disk is mostly optically thin to its 
own radiation, the midplane total emission dominates over 
the surface one (although the surface dust plays a crucial 
role for the heating of the disk). As a consequence, from our 
analysis we can extract information only on the dust opac- 
ity of the midplane, and for the rest of the paper when we 
consider the dust opacity we will mean only the midplane 
component. 

For the dust grains in the disk surface we assumed the 
same chemical composition and shape as for the midplane 
grains, whereas for their size distribution we assumed for 
a m i n and a max values around 0.1 pin (the sub-mm SED is 
insensitive to these values as long as a max in the surface is 
much lower than a max in the midplane, as expected if dust 
settling is important). 

Once the chemical composition and shape of the dust 
grains in the midplane are set, the dust opacity law depends 
on a m i n , a max and q. Considering a value of 0.1 /jm for a min 
(the dependence of the millimeter dust opacity from this 
parameter turns out to be very weak), the models param- 
eters for our analysis are the stellar parameters (L+, T e ff, 
M*), that have been set and listed in Table [2 plus the disk 
ones (Sdust.i; P, -Rout, q, a ma x), or equivalently (M dust , p, 
q, flmax, -Rout) where Md us t is the mass of dust in the disk. 

For reasons described in Section [5] it is convenient to 
approximate the millimeter dust opacity law discussed so 
far in terms of a power-law at millimiter wavelengths k = 
^imm(A/l mm) - ' 3 , with Ki mm m units of cm 2 per gram of 
dust. At a fixed value of q the relations between Ki mm , (3 
and a max can be determined. These are shown in Figure [3] 
for q = 2.5,3,3.5,4. 

4. Results 

The sub-mm/mm SED fits for all the objects in our sample 
are reported in Figure^B As discussed in Testi et al. {2003) 



We did not include existing NIR and mid-IR data in our 
analysis because at these shorter wavelengths the disk emission 



and Natta et al. (!2004[) . the only quantities that can be con- 
strained by the sub-mm/mm SED are the millimeter dust 
opacity spectral index (3 (in this paper calculated between 
1 and 3 mm) and the Mdust x Kimm producfl From the 
SED alone it is impossible to constrain the value of the 
surface density exponent p, since it is generally possible to 
fit the SED data with either very flat (p = 0.5) or very 
steep (p — 2) surface density profiles. However, all avail- 
able high angular resolution observations performed so far 
suggest that p < 1.5 (e.g. Dutrey et al. 119961 Wilner et 
al. I200D1 Kitamura et al. HE M Tes ti et al. I21KJ31 Andrew s 
& Williams [20071 Isella et al. [20091 Andrews et al. [20091) . 

The degeneracy of the SED on p and i? out has an im- 
pact onto the derived uncertainties for (3 and Md us t x «imm- 
Considering both this degeneracy and uncertainties of the 
observational data, the total absolute uncertainties are ap- 
proximately 0.2 — 0.3 for (3 and a factor of w 3 and w 4 
for Mdust x ^lmm for the spatially resolved and unmapped 
disks respectively. These uncertainties are mainly due to 
the adopted ranges for i? ut (listed in Table Q| and p be- 
tween 0.5 and 1.5, whose values determine the impact of 
the optically thick regions of the disk to the total emission, 
as explained in Section O 

The estimates for (3 and M dus t X ^imm obtained by our 
analysis are reported in Table 

4.1. Spectral slopes and dust opacity index 

Before analysing the values of /3, which provide information 
on the level of grain growth in the outer disk regions, it is 
important to check whether our results could be affected 
by detection biases: indeed, for a given flux at ~ 1mm, a 
disk with large dust grains in the outer parts will produce a 
stronger emission at ~ 3 mm than a disk with a population 
of smaller-sized grains and be easier to detect. In Figure [5] 
the grey region of the plot indicates the area in the (Fi mm , 
«i_3 mm ) plane in which the observations carried out so far 
are not sensitive, due to the sensitivity limits obtained at 1 
and 3 mm. While we could investigate the dust properties 
for 4 disks with 15 mJy < Fi mm < 30 mJy, the plot shows 
that a disk with a 1 mm flux in this same range, but with a 
spectral index ai_3 mm > 3.2, consistent with a disk model 
with ISM-like dust and an outer radius of ~ 100 AU, would 
not be detected at 3 mm. Nevertheless, we want to point 
out that only 2 of our 16 targets with Fi mm < 100 mJy 
were not detected at 3mm, and both have very weak -Fi.3 mm 
fluxes (upper limits of 27 mJy and 21 mJy for DP Tau and 
GK Tau respectively, Andrews & Williams 12005 j) . therefore 
the lack of disks with a ~ 3.2 at Fi mm > 30 mJy does not 
seem to be due to our sensitivity limit but is most likely 
real. Furthermore, as evident from the same plot, the 6 
faintest disks, with Fi mm < 60 mJy, show a spectral index 
between 1 and 3 mm that is lower on average (2.08 ± 0.08) 
than for the brighter disks (2.40 ± 0.06). This is the reason 
of the discrepancy of the average value of the dust opacity 



is optically thick and thus depends on the properties of the dust 
grains located in the surface layers of the inner disk regions. The 
aim of our analysis is instead to probe the optically thin disk 
emission in order to investigate the dust properties in the disk 
midplane. 

7 Note that, according to the adopted disk model, after fix- 
ing the stellar parameters, and these two quantities from the 
observed SED, the midplane temperature is well constrained. 
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Object 


iCt (K02) iCt (AW07) 


R b g (109) 


Adopted 7?out— interval 




(AU) (AU) 


(AU) 


(AU) 



AA Tau 




400l 6 7 °° 




200 


-400 


CI Tau 




225 ± 50 




150 


- 250 


CW Tau 








100 


-300 


CX Tau 








100 


-300 


CY Tau 


211 ± 28 




230 


230 


-330 


DE Tau 








100 


-300 


DL Tau 




175±§g 




150 


- 250 


DM Tau 


220 ± 47 


1501?™ 


160 


160 


- 260 


DN Tau 


147 ± 469 


iooir 


125 


125 


- 225 


DO Tau 


98 ±52 






100 


- 200 


DR Tau 


193 ± 141 


100+" 5 


90 


90 


- 190 


DS Tau 








100 


-300 


FM Tau 








100 


-300 


FZ Tau 








100 


-300 


GM Aur 




150 ± 25 


270 


270 


-370 


GO Tau 




350l«° 


160 


160 


- 260 


HO Tau 








100 


-300 


IQ Tau 


329l 4 7 °, 3 






250 


-400 


RY Tau 


81 ± 128 


150 ± 25 


115 


115 


- 215 


SU Aur 








100 


-300 


UZ Tau E 






160 


160 


- 260 



a) -Rout listed here results from image fitting using a truncated power-law for the surface density profile. K02: Kitamura et 
al. (|2002[l : AW07: Andrews & Williams (|2007|l . b) -R g is the radius within which the 95% of the source flux is observed, after 
fitting the visibility of the source with a gaussian; we used this value, when available, as a lower limit for our adopted 
-Rout-interval. 109: Isella et al. (|20"09j) . 



spectral index < (3 > between unmapped and resolved disks 
(see discussion below). 

In Figure [3] we plot Ki mm and (3 as functions of a max 
for different values of q. One interesting result is that the 
values of (3 obtained for the objects in our sample by the 
SED-fitting procedure are all less than or equal to 1, which 
can be explained only if dust grains with sizes larger than 1 
mm are present in the disk. For the 9 sources with /? < 0.5, 
the largest dust grains must have grown to sizes > 1 cm 
according to our dust model. If we let q assume values lower 
than 2.5, thus increasing the fractional number of the larger 
grains, the maximum grain size decreases, but a max >l mm 
for (3 < 1 is still needed. 

The estimates of the maximum grain sizes a max depend 
on the dust model that one adopts. If we consider for ex- 
ample the dust model adopted in Natta et al. (|2004|) (5% 
olivine, 15% organic materials, 35% water ice and 50% vac- 
uum) , the estimates for a max are greater by about an order 
of magnitude than the ones obtained with our dust model. 
In general, no realistic dust grain models result in (3 < 1 
for a m ax < 1 mm. 

Here it is worthwhile to remember that, whereas for the 
13 resolved disks in our sample (3 is well determined, for the 



8 disks for which we do not have information on their spatial 
extension, the optically thick parts of a possible compact 
disk could significantly contribute to the emission. If that is 
the case, no conclusions on the dust grain properties could 
be drawn. In particular, for the 6 unmapped disks with 
a ps 1.9-2.2 (CX Tau, DS Tau, FM Tau, FZ Tau, HO Tau, 
SU Aur) a very compact disk with i? ou t ~ 20 AU, a relative 
high inclination i pa 70°, and with ISM- like dust grains 
would fit the sub-mm/mm SED. For the other 2 unmapped 
sources, DE Tau (a ps 2.4) and CW Tau (a pa 2.5), ISM- 
like dust grains are consistent with the sub-mm/mm data 
using a disk model with i pa 70° and R out ps 30, 40 AU 
respectively. For the faintest disks, the outer radius has to 
be small because if the disks were larger they would also 
need to be more massive in order to keep being optically 
thick, thus producing an a value close to 2 with all the 
possible dust properties. In that case, they would emit more 
sub-mm/mm- wave radiation than observed. Since, for an 
optically thick disk, decreasing the inclination angle has 
the effect of making the SED steeper (because the central, 
hotter parts are less obscured), for disks with i < 70° the 
outer radii would need to be even smaller to fit the SED. 
Since these disk outer radii are lower by factors of 2 — 4 
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Fig. 3. Millimeter dust opacity for the adopted dust model. Left top panel: dust opacity at 1 mm as function of a max for 
a grain size distribution n(a) oc a~ q between a mul = 0.1 /im and a max . Different curves are for different values of q, as 
labelled. The dust grains adopted here are porous composite spheres made of astronomical silicates (« 10% in volume), 
carbonaceous materials (« 20%) and water ice (« 30%; see text for the references to the optical constants). Left bottom 
panel: (3 between 1 and 3 mm as a function of a max for the same grain distributions. Right panel: dust opacity at 1 mm 
as function of (3 between 1 and 3 mm for the same grain distributions; different iso-a max curves for a max = 0.1, 1, 10 cm 
are shown. In the plots the unit of measure for the dust opacity is cm 2 per gram of dust. 
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Fig. 5. Spectral index a between 1 and 3 mm plotted 
against the flux at 1 mm from the model best fits shown 
in Figure [4] The points plotted here are listed in the first 
two columns of Table [5] Open blue points represent the 
spatially resolved disks, whereas the filled red points are 
for the unmapped ones. A value of a consistent with a disk 
model with ISM-like dust and an outer radius of ~ 100 AU 
is indicated as a dashed horizontal line. The grey region 
shows the area in the plot that is unexplored because of 
the 3cr sensitivity limits of our new and past millimeter ob- 
servations (here considered as 15 mJy and 1.0 mJy at 1 and 
3 mm respectively). 



than the smallest disks that have been spatially resolved 
so far, and since for all the 13 resolved disks the values 
of (3 have been found always much lower than in the ISM 
(Asm ~ 1.7, see Figure [6]), thus showing clear evidence 
of grain growth, we consider very unlikely that many of 



the unmmapped sources are compact, optically-thick disks 
with unprocessed ISM grains. For AA Tau, which is the 
faintest of the mapped sources (Fi mm » 100 mJy) and has 
a low spectral index of a ~ 2.0 typical of fainter sources, 
the maps indicate a disk radius of 200 — 400 AU ruling out 
the possibility of optically thick disk regions dominating 
the total emission, and requiring a low f3 of 0.3 to explain 
the shallow slope. 

Also, only for the disks with (3 > 0.8, values of q > 3.5 
are consistent with the data (see right panel in Figure[3]), in 
which case a max > 10 cm. This is in good agreement with 
the following consideration using the approximated asyn- 
thotic formula (3 s» (q — 3)Prj, derived by Draine ([2006), 
where /3rj is the opacity spectral index of the solid material 
in the Rayleigh limit, and valid for a population of grains 
with 3 < q < 4 and a min <C a c <C a m ax, with a c being a 
function of the real and imaginary parts of the complex di- 
electric function and wavelength (at millimeter wavelengths 
a c ~ 1 mm). Since, at a fixed q, the function /3(a max ) de- 
creases when a max > 1 mm (see Figure [3]) , it follows that 
the approximated asynthotic value (q — 3)/3rj gives a lower 
limit for (3 at a finite a max . Adopting = 1.7 from the 
small dust grains in the ISM, for q>3.5 we derive /3>0.85, 
i.e. dust grains populations with q > 3.5 cannot explain j3- 
values lower than w O.8{0. 

In Figure [S] we report the histogram of the f3- values 
derived by our analysis. The average value in our sample 
considering both unmapped and resolved disks, < (3 >= 
0.6 ± 0.06, is not consistent with dust populations with q- 
values of 3.5, like for the ISM, or greater. Values of q lower 
than 3.5 have been obtained by, e.g., Tanaka et al. (2005) 
from simulations of grain growth in stratified protoplane- 
tary disks which give vertically integrated size distributions 



Note that /3-values greater than w 0.85 can be explained by 
both q < 3.5 and q > 3.5. 
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Fig. 4. Fits of the sub-mm/mm SEDs with the two-layer flared disk models (solid lines). Note that even if the data at 7 
mm, available only for 8 disks in our sample, are present in the plots, they have not been used in the fitting procedure 
(see Section [2751) . The errorbars in the plots take into account an uncertainty of 10% on the absolute flux scale, typical 
for flux measurements at these wavelengths. The circled data come from the new PdBI observations. The fitting values 
of the spectral index a and of the dust opacity spectral index /3 between 1 and 3 mm are indicated in the bottom left 
corner of each plot. The value of the surface density power- law index p adopted in these fits is 1, and the adopted value 
for the outer disk radius is listed in Table [5] 



with q w 3, plus a population of much larger bodies, to 
which however our observations would not be sensitive. 

These considerations clearly suggest that the process of 
dust grain growth to mm-sizes has played a fundamental 
role in most, if not all, of the Taurus- Auriga protoplane- 
tary disks whose sub-mm/mm SED has been sufficiently 
investigated. 

From Figure[6]it is also evident that the two sub-samples 
made of unmapped and resolved disks separately have dif- 



ferent averaged (3 values, with the /3 values being on av- 
erage lower for the unmapped disks than for the resolved 
ones (< (3 >= 0.41 ± 0.10, 0.68 ± 0.06 for the unmapped 
and resolved disks respectively) . Since the unmapped sam- 
ple nearly coincides with the fainter disks with Fi mm < 100 
mjy, and since (3 is correlated with the spectral index a, 
this is due to the trend already presented in Figure [5] be- 
tween -Fimm and ai_3 mm . This trend agrees with simple 
expectations for grain growth: as the maximum grain size 
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Table 5. Disk properties. 
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Object name 3 R^ ut _FT mm a^™ /3i_ 3mm M dust x Ki mm M^ t s M*~^ 

(AU) (mJy) (M x cmV 1 ) (M ) (M ) (M ) 



AA Tau 


300 


108 


2.0 


0.3 


2.8 


10" 


-4 


1.5 


10" 


-4 


1.5 


10" 


-3 








CI Tau 


200 


314 


2.5 


0.8 


7.9 


10" 


-4 


1.1 


10" 


-4 


1.2 


10" 


-4 


5.0 


10" 


-4 


CW Tau 


200 


129 


2.5 


0.8 


3.2 


10" 


-4 


4.2 


10" 


-5 


4.4 


10" 


-5 


1.9 


10" 


-4 


CX Tau 


200 


19 


2.3 


0.5 


5.7 


10" 


-5 


1.2 


10" 


-5 


1.9 


10" 


-5 








CY Tau 


280 


168 


2.2 


0.5 


6.0 


10" 


-4 


1.3 


10" 


-4 


2.0 


10" 


-4 








DE Tau 


200 


69 


2.5 


0.7 


1.6 


10" 


-4 


2.4 


10" 


-5 


2.8 


10" 


-5 








DL Tau 


200 


313 


2.4 


0.7 


8.4 


10" 


-4 


1.3 


10" 


-4 


1.5 


10" 


-4 








DM Tau 


210 


209 


2.5 


1.0 


9.2 


10" 


-4 


1.0 


10" 


-4 


1.0 


10" 


-4 


1.6 


10" 


-4 


DN Tau 


175 


153 


2.3 


0.5 


3.5 


10" 


-4 


7.7 


10" 


-5 


1.2 


10" 


-4 








DO Tau 


150 


220 


2.2 


0.4 


4.9 


10" 


-4 


1.4 


10" 


-4 


3.1 


10" 


-4 








DR Tau 


140 


298 


2.7 


1.0 


6.6 


10" 


-4 


7.4 


10" 


-5 


7.2 


10" 


-5 


1.2 


10" 


-4 


DS Tau 


200 


28 


2.0 


0.2 


6.2 


10" 


-5 


4.0 


10" 


-5 


1.1 


10" 


-3 








FM Tau 


200 


29 


1.9 


0.2 


6.5 


10" 


-5 


6.0 


10" 


-5 


1.7 


10" 


-3 








FZ Tau 


200 


27 


2.1 


0.2 


5.0 


10" 


-5 


3.2 


10" 


-5 


9.0 


10" 


-4 








GM Aur 


320 


423 


2.6 


1.0 


1.4 


10" 


-3 


1.6 


10" 


-4 


1.5 


10" 


-4 


2.5 


10" 


-4 


GO Tau 


210 


151 


2.6 


1.0 


6.3 


10" 


-4 


7.1 


10" 


-5 


6.9 


10" 


-5 


1.1 


10" 


-4 


HO Tau 


200 


36 


2.1 


0.5 


1.7 


10" 


-4 


3.7 


10" 


-5 


5.6 


10" 


-5 








IQ Tau 


325 


118 


2.2 


0.5 


3.4 


10" 


-4 


7.4 


10" 


-5 


1.1 


10" 


-4 








RY Tau 


165 


383 


2.3 


0.5 


5.1 


10" 


-4 


1.1 


10" 


-4 


1.7 


10" 


-4 








SU Aur 


200 


50 


2.1 


0.2 


5.0 


10" 


-5 


3.2 


10" 


-5 


9.0 


10" 


-4 








UZ Tau E 


210 


333 


2.5 


0.7 


6.1 


10" 


-4 


9.6 


10" 


-5 


1.1 


10" 


-4 









a) Underlined objects are those that have not been mapped to date through high-angular resolution imaging, b) Central value of 
the adopted R ou t -interval reported in Table [4] This is the value used to extract the other quantities listed in the Table, c) Source 
flux density at 1 mm from the best fit disk model described in Section [3] 



a max increases, both /3 and K lmm decrease, so disks with the 
same dust mass should become fainter as (3 gets smaller. 
Alternatively, the Fi m m-Qi-3mm trend could arise if the 
fainter disks are smaller, and then more optically thick, 
than the brighter, resolved ones, while their "true" (3 would 
be nearly the same. 

In Figure [7] we show the (3 vs stellar age plot. No ev- 
idence of an evolutionary trend is present over the large 
range of stellar ages spanned by our sample. As already 
discussed the determination of the individual ages of PMS 
stars is very uncertain. Hartmann (|200ip has shown that 
the large age spread deduced for T Tauri stars in Taurus 
may just be the result of uncertainties in the measurements 
towards individual members of a population with a very 
narrow age spread. Our most solid conclusion, not affected 
by the age determination uncertainties, is that most (if not 
all) of the Class II YSO in our sample show evidence of 
evolved grains in the disk outer regions. 

4.2. Grain growth 

These results confirm two statements that have been ad- 
dressed in the last years. The first one is that the dust 



coagulation from ISM dust to mm/cm-sized grains is a fast 
process, that appears to occur before a YSO enters in the 
Class II evolutionary stage, i.e. when a dense envelope still 
surrounds the protostar+disk system. This is also consis- 
tent with the recent results by Kwon et al. p00 9i) who have 
found evidence of grain growth in 3 Class YSOs in Perseus 
and Cepheus through CARMA observations at 1.3 and 2.7 
mm. 

The second statement is that the mm/cm-sized dust 
grains appear to stay longer in the outer disk if com- 
pared with their expected inward drift timescale due to 
the interaction with the gas component. In a circum- 
stellar disk the gas experiences a pressure gradient force 
which cancels a part of the central star's gravity and in- 
duces a rotation slower than the Keplerian velocity. Dust 
grains, whose orbits are hardly affected by gas pressure, 
are assumed to be rotating with Keplerian velocity, and 
thus experience a headwind. They lose angular momentum 
and spiral inward to the central star (Adachi et al. [1976, 
Weidenschilling fl977p . Takeuchi & Lin (2005) proposed as 
possible solutions to the drift problem that either the grain 
growth process from mm-sized to cm-sized grains is very 
long in the outer regions of the disk (timescales greater than 
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Fig. 6. Distribution of the dust opacity spectral index (3 for 
the disks in our sample. In blue the values for the spatially 
resolved disks are indicated, whereas the red is for the un- 
mapped ones. The average [3 value for all the sources in 
our sample (black), for the unmapped sources only (red), 
for the resolved sources only (blue), and the value of 1.7 for 
the ISM dust are indicated as dashed vertical lines. 



2.0 
1.5 
^ 1.0 
0.5 
0.0 



• unmapped TTS 
o resolved TTS 



0.1 



1.0 



10.0 



age 



(Myr) 



Fig. 7. Beta versus age: relationship between the dust opac- 
ity spectral index (3 and the estimated stellar age obtained 
as in Section [2.31 Open blue points represent the spatially 
resolved disks, whereas the filled red points are for the un- 
mapped ones. The (3 value for the ISM dust is indicated as 
a dashed horizontal line. 



1 Myr, e.g. given by grains' sticking probability as low as 
0.1) or that grains have already grown to 10 m or larger and 
they do not migrate rapidly anymore, and mm-sized parti- 
cles are continuosly replenished through collisions of these 
large bodies. However our analysis suggests that, unless the 
grain size distribution index is very low, i.e. q < 2.5, disks 
with (3 < 0.5 show evidence of grains as large as ~ 1 cm in 
the outer regions. Also, as already pointed out in Natta et 
al. (|2004p and Brauer et al. (20(37|), a noteworthy amount of 
larger bodies, which would not contribute to the disk emis- 
sion because of their negligibly small opacity coefficient, 
could significantly increase the disk masses. Since the typi- 



cal orbital decay time for a cm-sized grain with a density of 
0.1 g cm -3 at 50 AU is ps 10 5 yr (10 4 yr for a more compact 
grain with a 1 g cm -3 density), to explain the presence of 
cm-sized grains in the older disks in our sample (with ages 
of few Myrs and up to ~ 10 7 yr) the reservoir of larger bod- 
ies would have to contain at least 100 times more mass than 
the observed dust mass. The total disk mass would then be 
comparable to or even exceeding the estimated stellar mass 
(using a standard gas-to-dust mass ratio of 100, see Section 
EM- 

To solve the drift problem for very low and very high 
disk masses Brauer et al. (|2007p proposed a significant re- 
duction of the drag experienced by the dust due to a reduc- 
tion of the gas-to-dust ratio in the disk from the canonical 
value of 100, and to possible collective effects of dust gath- 
ering in a thin midplane layer in the case of very-low tur- 
bulent disks. For example, using a low value for the turbu- 
lence parameter a — 10~ 6 and a gas-to-dust ratio of 5, they 
find that cm-sized dust grains remain in the outer parts of 
the disk for more than 2 Myr if disk masses < 0.05 M* or 
> 0.2 are considered. In the range Mdisk ~ 0.05—0.2 
the collective effects of dust are not efficient enough to re- 
tain cm-sized grains in the disk outer regions. If the disk is 
more turbulent, then more gas has to be removed in order 
to keep the particles in the outer parts for a longer period 
of time. To justify such small gas-to-dust ratios tha authors 
suggest the process of photoevaporation as the main driver 
of gas removal. However, in this case, since photoevapo- 
ration is expected to manifest only in the later stages of 
disk evolution, then we would expect a strong evolutionary 
trend of (3, which is instead not found by our analysis. 

Another possibility is the presence in the disk of local 
gas pressure maxima that would block the inward drift of 
the dust grains. As shown by the dynamics calculation of 
Barge & Sommeria (|1995[) and Klahr & Henning (|1997p this 
mechanism may provide a solution also to the relative veloc- 
ity fragmentation barrier that prevents meter-sized bodies 
to further grow in size in the inner regions of the disk. 

4.3. Dust mass 

From the sub-mm/mm SED fitting procedure we derived 
the quantity Afdust x Ki mm . To obtain the dust mass in 
the disks an estimate of /ti mm is needed. The choice of the 
dust opacity at 1 mm creates the largest uncertainty in the 
derivation of (dust) disk masses from millimeter observa- 
tions. Estimates of Ki mm from different dust models in the 
literature span a range of a factor greater than 10. The 
most commonly used value in the literature is Ki mm ~ 3 
cm 2 per gram of dust (Beckwith et al. 1990) that, for sim- 
plicity, is supposed to be constant, i.e. independent of the 
dust grain properties, for all the disks. As noted in Natta 
et al. (|2004p this is not realistic as grain growth affects 
both the frequency dependence of the dust opacity and its 
absolute value (e.g. Miyake & Nakagawa ll993p . 

In this paper we use the dust model described in Section 
13. II to get an estimate of Ki mm that takes into account the 
dust physical properties, expecially grain growth. As shown 
in the left top panel of Figure^ at a fixed value of the grain 
size distribution index q, the value of Ki mm can be inferred 
from a max which in turn can be derived from (3 through the 
relation plotted in the left bottom panel of the same figure. 
The relation between /ti mm and (3 is shown in the right 
panel. Therefore, from an estimate of (3 and Afdust x Kimm 
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obtained by the sub- mm/mm SED fitting procedure we can 
get an estimate for Ki mm and then the dust mass Afd us t for 
our dust model, after fixing the q parameter. In Table [5] we 
report the values of Md us t as obtained with this method for 
q = 2.5, 3.0, 3.5. It is important to bear in mind that, as 
for the estimate of the maximum grain size, also for Ki mm , 
and thus for the dust mass, the estimates depend on the 
specific dust model that one adopts (see discussion in Natta 
et al. 12001]) . However, following Isella et al. (|2009)) . we take 
into account the variation of the dust opacity with (3 using 
a well defined dust model to infer dust masses, and this 
approach is more accurate than considering a fixed rei mm 
as often done in previous studies. 

In general, for our specific dust model, when (3 < 1 
then M|^- 5 < < Mj^' 5 , the amount of variation 

of M| ust for the different q values being dependent on (3. 
In particular, when (3 ~ 0.8 — 1 M%~ st ' and M d ^ t axe 

nearly coincident and M%~ st ' is greater than M%~ st by 
factors of ~ 2 — 4. For (3 < 0.7 a grain size distribution 
with 5 = 3.5 cannot explain the observational data and the 
discrepancy between M^ t 5 and increases up to a 

factor of ~ 15 when (3 ~ 0.2. This is due to the increase of 
the distance between the Ki mm ((3) curves at fixed q- values, 
when (3 decreases, as shown in the right panel of Figure [3] 

As already mentioned in Section 14.11 from the value of 
(3 it is impossible to disentangle a max and q. Nevertheless, 
only the disks showing (3 > 0.8 can be explained by a grain 
size distribution index q = 3.5. For this reason, for the rest 
of our analysis we will consider separately only the cases 
in which the dust grain size distribution is characterized 
by q = 2.5 or 3. The value of q in protoplanetary disks is 
unknown and it is not observationally constrainable, and 
so it is in principle possible that it can vary substantially 
from one disk to the other. For our discussion we will use 
the assumption that q is the same for all the disks. 

Once the dust mass is estimated, this is commonly con- 
verted into an estimate of the total disk mass (dust+gas) 
using the standard ISM gas-to-dust ratio of 100. However, 
if the dust and gas components have different evolutionary 
timescales, this ratio is a function of time and the stan- 
dard ISM value would be inappropriate for estimating the 
disk mass at the present time. Physical processes that are 
likely to alter this ratio during the evolution of a proto- 
planetary disk are gas photo-evaporation, which would lead 
to a decrease of the gas-to-dust ratio (unless a significant 
amount of dust is dragged away by the gas flow), and dust 
inward migration, which instead would increase it. Natta et 
al. ((2004) suggested that the standard ISM ratio provides a 
correct estimate of the "original" disk mass, when the disk 
composition reflected that of the parent cloud. However, 
this is true only if a negligible amount of dust has migrated 
with respect to the gas and accreted onto the central star, 
otherwise the ISM ratio would provide only a lower limit 
to the original disk mass. Because of our ignorance on the 
gas-to-dust ratio evolution in this paper we will not convert 
the dust mass into a total disk mass. Also, it is important 
to note that since millimeter observations are insensitive 
to the emission of bodies much larger than 1 cm, including 
planetesimals and planets, the dust mass we are referring to 
in this paper does not include the contribution from these 
larger objects. The dust mass is intended to be the mass in 
"small grains" only, and is a lower limit to the total mass 
in dust. 



Figure [5] shows the (3 vs Mj us t plots for q equal to 2.5 
and 3. The plot with q = 2.5 does not show any signifi- 
cant relation. In the q = 3 case an anticorrelation between 
these two quantities may be present, with the more massive 
dust disks being associated to the ones showing today the 
largest grains in their outer parts, with a max up to ~ 10 
cm. This possible anticorrelation is clearly an effect of the 
choice of a physical model for the dust that takes into ac- 
count the variation of the dust opacity normalization factor 
with the size-distribution of the dust grains (Section [BTTJ) ■ 
If we had made the unrealistic assumption of a constant 
value for the dust opacity at, say, 1 mm, we would have 
found an opposite result, i.e. disks with less mass in dust 
would have shown on average larger grains, as can be de- 
rived from Figure [5l since a is approximately proportional 
to (3 and -Fi mm would be approximately proportional to 
the dust mass, in the constant Ki mm case. More realisti- 
cally, lower values of (3 imply low absolute values of the 
dust millimeter opacity (see right panel of Figure [3]) , so 
that, for a fixed observed flux and disk temperature, the 
corresponding dust mass must be larger. The reason why 
the plot with q = 2.5 does not show signs of anticorrela- 
tion is that at a fixed (3 a lower q value is associated to a 
smaller dust grain maximum size. The f3 range spanned by 
the disks in our sample translates into a range in Ki mm that 
is narrower for q = 2.5 than it is for q = 3 (see right panel 
of Figure [3]), thus leading to smaller correction factors with 
respect to the constant Ki mm case. Note that among the 5 
disks showing the smallest (3 (« 0.2 — 0.3) and the highest 
dust masses (about 10 Mq for q = 3) only one (AA Tau) 
has been mapped and spatially resolved. High angular res- 
olution and high sensitivity imaging of these disks together 
with other ones are needed to test the robustness of the 
(3 — A/disk, init anticorrelation in the q = 3 case. 

Figure [5] and [TU] plot the disk mass in dust against the 
estimated mass of the central star (dust mass obtained 
with q-values of 2.5 and 3) and mass accretion rate (dust 
mass with q = 3 and stellar ages smaller and greater 
than 2 Myr) respectively. In the 0.3 — 2 Mq range in 
stellar mass the dust mass does not show any trend (for 
q = 2.5: 3 • 10~ 5 < M dust < 2 • 10- 4 M Q ; for q = 2.5: 
4 • 10~ 5 < M disk ,mit < 2 • 1O~ 3 M ). The two YSOs with 
M* < 0.3 Mq have lower dust masses than all the other 
YSOs for both q — 2.5 and 3. However, the range in stellar 
mass spanned by our sample is rather limited, and partic- 
ularly biased toward YSOs with late-K and early-M spec- 
tral types. From the upper limits obtained at 1.3 mm by 
Schaefer et al. ([2009) for 14 T-Tauri stars with spectral 
types later than ss M2 (i*i. 3mm < 5 — 20 mJy) , and from the 
detection of very low disk masses around sub-stellar objects 
in Taurus- Auriga (Klein et al. [20031 Scholz et al. 12006ft , it is 
likely that a correlation between the disk mass in dust and 
the mass of the central star exists. High sensitivity millime- 
ter observations are needed to constrain the disk masses in 
dust around very low mass PMS stars. In Figure [TUl a trend 
of higher mass accretion rates for the more massive disks 
in dust may be present in the data but still the scatter 
is very large (a similar result is obtained for the q — 2.5 
case). Other relations between the dust properties {(3, dust 
mass) and the stellar ones (e.g. stellar luminosity, effective 
temperature) have been investigated but no significant cor- 
relations have been found. 

Finally, we have compared the values of (3 obtained in 
this paper with the ones derived in Natta et al. (2004) for 
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Fig. 8. Dust opacity spectral index (3 and maximum grain size a max plotted against the disk mass in dust for the sources 
in our sample. In the left and right panels the dust mass has been obtained adopting a value of 2.5 and 3 respectively for 
the power index q of the dust grains size distribution. Open blue points represent the spatially resolved disks, whereas 
the filled red points are for the unmapped ones. The errorbars in the dust mass estimates do not take into account our 
uncertainty on the real value of q. 
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Fig. 9. Disk mass in dust plotted against the stellar mass for the sources in our sample. In the left and right panels 
the dust mass has been obtained adopting a value of 2.5 and 3 respectively for the power index q of the dust grains 
size distribution. The stellar mass has been estimated as described in Section I2.3[ and the errorbars take into account 
the range of values given by different PMS stars evolutionary models. Open blue points represent the spatially resolved 
disks, whereas the filled red points are for the unmapped ones. The errorbars in the dust mass estimates do not take into 
account our uncertainty on the real value of q. 



a sample of 6 intermediate mass Herbig Ae stars and 3 
lower mass T Tauri stars. From the two-sided Kolmogorov- 
Smirnov test, the probability that the values from the two 
samples come from the same distribution is about 88%, in- 
dicating that the disks in the two samples have similar dust 
properties in terms of grain growth (in the Natta et al. 2004 
sample only one object, the Herbig Ae star HD150193, has 
a relatively high (3 = 1.6±0.2, compatible with non-evolved 
ISM- like dust grains). 



5. Summary 

We have analysed the sub-mm/mm SED out to 3mm of a 
sample of 21 Class II T- Tauri stars in the Taurus- Auriga 
star forming region. This sample comprises about the 60% 
of the "isolated" class II YSOs with an estimated stellar 
mass greater than w 0.4 M©. It also comprises approxi- 
mately 1/3 of the isolated class II wich have been observed 
at 0.85mm to have a flux less than 100 mJy, for which we 
provide first 3mm detections thanks to sensitive PdBI ob- 
servations. Our main findings are summarized below: 

1 . For all the sources in our sample the millimeter opacity 
spectral index (3 is significantly lower than the value 
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Fig. 10. Disk mass in dust plotted against the mass accretion rate for the sources in our sample. In the left and right 
panels we included the objects with estimated ages of less and more than 2 Myr respectively. The dust mass has been 
obtained adopting a value of 3 for the power index q of the dust grains size distribution. The stellar mass accretion 
rate has been estimated following the method described in Section [2~4l Open blue points represent the spatially resolved 
disks, whereas the filled red points are for the unmapped ones. The errorbars in the mass accretion rate take into account 
the typical intrinsic time-variability of the accretion rates. The errorbars in the dust mass estimates do not include our 
uncertainty on the actual value of q. 



obtained for the dust in the ISM. For the 13 spatially 
resolved disks this is a clear evidence of the presence in 
the outer disk regions of dust grains as large as at least 1 
mm. For the 8 unmapped sources, assuming disk outer 
radii comparable to the resolved sources low (3 values 
are found as well, suggesting that even for these fainter 
sources 1 mm-sized grains are in the disk outer regions. 
This confirms, over a larger sample and less massive 
disks, the results obtained by past observations on T- 
Tauri stars (e.g. Rodmann et al. 120061 and references 
therein) . 

2. No significant evidence of an evolutionary trend for the 
millimeter opacity spectral index j3 has been found. This 
indicates that the dust grain growth to ~ 1 mm-sizes is a 
very fast process in a protoplanetary disk, that appears 
to occur before a YSO enters in the Class II evolutionary 
stage. Also, the amount of these large grains in the disk 
outer regions does not appear to decline throughout all 
the Class II evolutionary stage. 

3. Only for 9 sources in our sample a power law index for 
the grain size distribution of the ISM, <7ism ~ 3.5, is 
consistent with the sub/mm-mm SED. Instead, most of 
the disks show evidence of a smaller q— value. 

4. We have not found any significant correlation between 
the dust properties, namely grain growth and dust mass 
in the disk outer regions, and the properties of the cen- 
tral star, including the mass accretion rate onto the stel- 
lar surface. However, our sample contains mostly YSOs 
with a narrow range in spectral types, namely between 
late-K and early-M, and it is possible that a correlation 
between the disk mass in dust and the stellar mass will 
become apparent when considering also disks around 
late-M spectral types. 

5. The 6 faintest sources in our sample, with F lmm < 50 
mjy, show a spectral index between 1 and 3 mm that is 
on average lower than the spectral index of the brighter 
sources. This may indicate either that for these fainter, 



yet unmapped disks the emission from the optically 
thick inner disk is much more significant than for the 
brighter sources or that the dust grains in the outer 
regions of these fainter disks are even larger than in 
the brighter YSOs, possibly up to 10 cm if q = 3. The 
latter hypothesis would imply on average a lower abso- 
lute value for the millimeter opacity coefficient for these 
fainter disks. Taking into account this effect through a 
physical model for the dust properties, we found, in the 
case of a value for the power index of the grain size 
distribution q = 3, an anticorrelation between (3 and 
the disk dust mass, that would indicate that larger dust 
grains are found in the outer parts of the more massive 
disks (in dust). 

In the next future, with the advent of the next sub- 
mm/mm facilities (both new interferometers like ALMA 
and technical improvements in the existent arrays, like for 
example the E-VLA) it will be possible to extend the in- 
vestigation of the dust properties to fainter circumstellar 
disks around the lowest mass PMS stars and brown dwarfs. 
Furthermore, thanks to the great enhancements in both 
high angular resolution and sensitivity that these facilities 
will provide, it will be possible to determine their spatial 
extension. 
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